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ECIFICATION 



METHOD FOR MEASURING DIGGING POSITION 




BACKGROUND OF THE INVENTION 
[Technical Field of the Invention] 

The present invention relates to a method for determining the digging 
position in a non-open-cut method of excavation and, more particularly, to a 
5 method which ensures accuracy in determining positions by lessening the 
influence of a noise magnetic field of frequency components close to that of the 
signal magnetic field to be measured. 

(Prior Art] 

A horizontal drilling method, which is one of the non-open-cut methods of 
10 this kind, uses a small-diameter pipe of 100 mm or less across for horizontally 
digging in the ground, and accordingly, such a precision position detennining 
apparatus as used in an ordinary small-diameter driving method of excavation 
cannot be placed near a drill. To solve this problem, it is customary in the art to 
adopt a method in which an AC magnetic field is generated by a coil mounted in 
15 the drill head and detected by an above-ground magnetic sensor like a coil to 
^ determine the current digging position. 

This method is simple and easy, but since the magnetic field by the coil is 
a dipole magnetic field, it rapidly attenuates with distance from the coil. Hence, 
this method has a defect of inability to achieve high-reliability determination of the 
20 digging position when a power line or similar magnetic noise source is present in 
the vicinity of the place where to perform the position determination. 

SUMMAR OF THE INVENTION 
An object of the present invention is to provide a position determination 
method that permits high-reliability determination of the digging position by 
25 detecting on the ground the AC signal magnetic field provided from a coil housed 
in the drill head even if a noise magnetic field is present which affects the position 
determination. 

To attain the above object, a digging position determining method 
according to the present invention for non-open-cut excavation, which senses an 
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AC magnetic field provided from a magnetic field source by a magnetic sensor 
provided on the ground and calculates the position of the magnetic field source 
from the magnitude and direction of the sensed magnetic field, said method having 
a construction characterized in that 

In case where in addition to a signal magnetic field generated by said 
magnetic field source, there exists a noise magnetic field generated by a nearby 
current, 

at lease one of the position of said magnetic field source, the tilt angle of 
said magnetic field source to the vertical direction and the azimuth angle of its 
axial direction of said magnetic field source in a horizontal plane is calculated, 
from a projective component of the magnetic field sensed by said magnetic sensor 
and projected on a plane or straight line orthogonal to a vector-valued direction of 
said noise magnetic field. 

That is, in 1999 year's investigation and research relating to useful 
utilization techniques of energy resources: entitled "Research for low-loss 
energizing techniques in establishment of advanced telecommunication 
network",(executed by Composite Development System for New Energy 
Industrial Technique), it is described that external noise magnetic fields, which 
affects the position determination in the non-open-cut method of excavation, is 
mostly generated by a current of some kind. In this case, the magnitude of the 
noise magnetic field varies irregularly with time, but its vector-valued direction is 
constant at each field sensing position. 

The present invention attains its object through the adoption of the 
following steps (A) and (B). 

(A) The direction of the noise magnetic field is detected, and a sensed 
magnetic field in which the noise magnetic field and a signal magnetic field are 
mixed is projected on a plane or straight line orthogonal to the direction of the 
noise magnetic field to obtain a projective component. 

Since the projective component is theoretically free from a component 
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derived from the noise magnetic field, at lease one of the position, azimuth angle 
and tilt angle of the magnetic field source is calculated so that the magnitude of the 
projective component (in the case of projection on the straight line) or it magnitude 
and direction (in the case of projection on the plane) is substantially equal to a 
5 theoretically calculated value of a corresponding quantity of a magnetic field 
generated from a magnetic field source, or ,has a minimum difference between the 
former and the latter. The sensed magnetic field is obtained by sensing magnetic 
fields at different positions whose number is determined by how many ones of the 
position, azimuth angle and tilt angle of the magnetic field source are unknowns 
1 0 and how such unknowns are calculated. 

(B) To obtain the noise magnetic field, 

a) When the number of noise magnetic fields is virtually one: 

• The noise magnetic field is sensed to obtain its direction essentially in the 
absence of a signal magnetic field. 
15 • When the noise magnetic field has frequency components at frequencies 

different from that of the signal magnetic field, the frequency components are 
measured to obtain the direction of the noise magnetic field 

b) When the number of magnetic fields is virtually two: 

The frequency components of a first one of the two noise magnetic fields, 
20 which are widely spaced from the frequency components of the second noise 
magnetic field and the signal magnetic field, are measured to obtain a vector 
direction of the first noise magnetic field, and the frequency components of the 
second noise magnetic field, which are widely spaced from the frequency 
components of the first noise magnetic field and the signal magnetic field, are 
25 measured to obtain a vector direction of the second noise magnetic field. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The features of the present invention will be clearly understood from the 
following description taken in conjunction with the accompanying drawings, in 
which: 
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Fig. 1 is a perspective view explanatory of the placement of a magnetic 
sensor in the present invention method; 

Fig. 2 is a vector diagram explanatory of the principle of measurement 
according to the present invention method in a case of one noise magnetic field 
source; 

Fig. 3 is a vector diagram explanatory of the principle of measurement 
according to the present invention method in a case of two noise magnetic field 
sources; 

Fig. 4 is a flowchart illustrating a procedure for obtaining a projective 
component according to the present invention method in a case of one noise 
magnetic field source; 

Fig. 5 is a flowchart illustrating a procedure for obtaining a projective 
component according to the present invention method in a case of two noise 
magnetic field sources; 

Fig. 6 is a flowchart illustrating a procedure for obtaining a signal magnetic 
field in a case of one noise magnetic field; 

Fig. 7 is a flowchart illustrating a procedure for obtaining a signal magnetic 
field in a case of two noise magnetic fields; 

Fig. 8 is a flowchart showing a measurement procedure in the present 
invention method; 

Fig. 9 is a flowchart showing a procedure for calculating the position of the 
signal magnetic field according to the present invention method in a case of one 
noise magnetic field when the number of unknowns and the number of equations 
are the same as to each other- 
Fig. 10 is a flowchart showing a procedure for calculating the position of 
the signal magnetic field source according to the present invention method in a 
case of one noise magnetic field source when the number of equations is larger 
than the number of unknowns; 

Fig. 11 is a flowchart showing a procedure for calculating the position of 




the signal magnetic field according to the present invention method in a case of 
two noise magnetic field sources when the number of unknowns and the number 
of equations are the same as to each other; 

Fig. 12 is a flowchart showing the procedure for calculating the position of 
5 the signal magnetic field source according to the present invention method in a 
case of two noise magnetic field sources when the number of equations is larger 
than the number of unknowns; 

Fig. 13 is a flowchart showing an example of a procedure for obtaining the 
signal magnetic field through the use of the magnitude of the signal magnetic field 
10 vector; 

Fig. 14 is a flowchart showing another example of a procedure for 
obtaining the signal magnetic field through the use of the magnitude of the signal 
magnetic field vector; 

Fig. 15 is a perspective view explanatory of the placement of a magnetic 
15 sensor for determining the direction of a noise magnetic field according to the 
present invention method when a digging head, which is a signal magnetic field 
source, is distant from the position of measurement; 

Fig. 16 is a flowchart showing a procedure for calculating the direction of 
the noise magnetic field according to the present invention method in a case where 
20 the digging head is distance from the position of measurement and the signal 
magnetic field source stops to generate the magnetic field; 

Fig. 17 is a flowchart showing a procedure for calculating the direction of 
the noise magnetic field according to the present invention method when the signal 
magnetic field and the noise magnetic field are mixed to each other; 
25 Fig. 18 is a signal frequency spectrum diagram explanatory of how the 

operation of selecting a frequency for maximizing a frequency spectrum is used in 
the process for calculating the direction of the noise magnetic field in the flowchart 
shown in Fig. 13; 

Fig. 19 is a flowchart explanatory of a first method by which the flow of 



candidate vector calculating process is used in the process for calculating the 
direction of the noise magnetic field in the flowchart shown in Fig. 17; 

Fig. 20 is a flowchart explanatory of a second method by which the flow of 
candidate vector calculating process is used in the process for calculating the 
5 direction of the noise magnetic field in the flowchart shown in Fig. 1 7; 

Fig. 21 is a flowchart explanatory of a method by which the flow of 
process for evaluating the candidate vector and obtaining the direction of the noise 
magnetic field is used in the process for calculating the direction of the noise 
magnetic field in the flowchart shown in Fig. 17; 
10 Fig. 22 is a signal frequency spectrum diagram explanatory of how the 

operation of selecting a frequency for maximizing a frequency spectrum is used in 
the process for calculating the direction of the noise magnetic field in the flowchart 
shown in Fig. 13. 

Fig. 23 is a signal waveform diagram explanatory of the operation of 
15 specifying a period in which only a noise magnetic field exists through utilization 
of a fact that the amplitude of the sensed magnetic field signal becomes small 
during the OFF period of the signal magnetic field that is turned OFF by a 
predetermined procedure, in the calculation of the direction of the noise magnetic 
field according to the present invention method; 
20 Fig. 24 is a signal waveform diagram showing instantaneous variations in 

the amplitude of the sensed magnetic field signal when the signal magnetic field is 
periodically turned OFF in the calculation of the direction of the noise magnetic 
field according to the present invention method; 

Fig. 25 is a signal waveform diagram showing instantaneous variations in 
25 the amplitude of the sensed magnetic field signal when the signal magnetic field is 
randomly turned OFF in the calculation of the direction of the noise magnetic field 
according to the present invention method; 

Fig. 26 is a flowchart showing a method for calculating the direction of the 
noise magnetic field according to the present invention in which the signal 
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magnetic field is turned OFF on the basis of a predetermined sequence, a period 
during which a particular magnetic field is OFF is specified by a sequence starting 
at the time when a correlation function between the sequence and the sensed 
magnetic field becomes maximum, and the direction of the sensed magnetic field 

5 in the specified period is regarded as the direction of the noise magnetic field; 

Fig. 27 is a flowchart showing a method for calculating the direction of the 
noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and the 
starting time of a sequence indicating the most likely ON/OFF state of the signal 

10 magnetic field is calculated from a plurality of times when a correlation function 
between the predetermined sequence and the sensed magnetic field becomes 
maximum; 

Fig. 28 is a flowchart showing the method for calculating he direction of 
the noise magnetic field according to the present invention in which the signal 
15 magnetic field is turned OFF on the basis of a predetermined sequence, the 
correlation function between the predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points at which the period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variance is 
20 minimum is regarded as the direction of the noise magnetic field; 

Fig. 29 is a flowchart showing another method for calculating the direction 
of the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, a 
correlation function between the predetermined sequence and the sensed magnetic 
25 field is calculated at each of a plurality of time points at which the period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variance is 
minimum is regarded as the direction of the noise magnetic field; 

Fig. 30 is a flowchart showing the method for calculating the direction of 
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the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and the 
starting time of a sequence indicating the most likely ON/OFF state of the signal 
magnetic field is calculated from a plurality of time points when the correlation 
5 function between the predetermined sequence and the sensed magnetic field 
becomes maximum; 

Fig. 31 is a flowchart showing still another method for calculating the 
direction of the noise magnetic field according to the present invention in which 
the signal magnetic field is turned OFF on the basis of a predetermined sequence, a 

10 correlation function between the predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points into which the period of the 
sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variance is 
minimum is regarded as the direction of the noise magnetic field; and 

15 Fig. 32 is a perspective view depicting an example of a magnetic field 

sensing frame for use in the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
As depicted in Fig. 1, in a case where the position of a digging head 2 
under the ground surface 1, which is a signal magnetic field source to be sensed, is 

20 determined using a magnetic sensor 4 placed on the ground surface 1 at a proper 
position when a power line or similar magnetic noise source 3, which generates a 
noise magnetic field, is placed near the digging position to be determined, there are 
present a signal magnetic field vector Hs provided from the digging head 2 and a 
noise magnetic field vector Hn from the magnetic noise source 3, such as power 

25 line. In this case, the magnetic sensor 4 senses a vector H m that is a combined 
version of the signal magnetic field vector Hs and the noise magnetic field vector 
H n . 

Now, let the noise magnetic field of a position vector r and at time t be 
identified as a vector H„(r, t). On the other hand, let the signal magnetic field 
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generated by magnetic filed generating means for position sensing be identified as 
a vector Hs(r~r c , 0 C , t). Here, the vector 9 C is an angle of orientation of the 
magnetic field generating means, which is defined by three angles of rotation in the 
coordinate system fixed to the ground that is the coordinate system fixed to the 
magnetic field generating means. 

Since the noise magnetic field and the signal magnetic field are both 
sensed simultaneously and since the noise magnetic field vector H^t, t) varies 
randomly with time, it is impossible to extract only the signal magnetic field vector 
H s (r-r c , 0 C , t) from the measured magnetic field vector H m (r-re, 0 C , t) unless 
statistical properties of the noise magnetic field are known and the noise magnetic 
field is signal-wise orthogonal to the signal magnetic field Even if the statistical 
properties for separating the noise magnetic field from the signal magnetic field are 
known prior to the determination of the position of the latter, the separation calls 
for a large amount of data, and hence the conventional scheme is not ever practical. 
According to the present invention, the direction of a vector e^r) of the noise 
magnetic field vector H n (r, t) is obtained by separated means, and in a coordinate 
system shown in Fig. 2, a component vector H m p (r-r c , G c , t), projected to a plane 
vertical to the direction of a vector e^r) of the measured magnetic field vector 
H m (r-r c , 0 C , t) as shown in Fig. 4 (S 1, S2, S3). 

H m P (r-r c , 0 C , /)=H m (r-r c , 0(H m (r-r c , & /> e n (r)> n (r> (1) 
This component does not contain the noise magnetic field for the reason given 
below. Since 

H m p (r-r c , 0 C , 0=U s (r- r c , 0 c , /)+H n (r, 0=H s (r-r c , 3, /)+|H n (r, 0|e n (r). (2) 
it follows that 

H m P (r-r c , 6>, /)= H s (r- r c , 0 Q9 0"(H 8 (r-r c , 0 C , /). e n (r)> n (r). (3) 
from which it can be seen that the projective component vector H m p (r-r c , 6 C , 0 does 
not contain the component of the noise magnetic field vector H„(r, t). 

However, the projective component vector H m p (r-r c , 9 C , t) has lost 
information of one axis by the projection on a plane vertical to the direction of the 
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vector e„(r). That is, since the same projective components are obtained 
irrespective of the magnitude of a component parallel to the vector en(r), two 
independent components are obtained. 

Although any method can be used to obtain the two independent 
5 components; it is possible to use such a method as described below. 

That one of coordinate axes of a measurement coordinate system C M 
(which will be described later on), which is not parallel to the direction e„(r) of the 
noise magnetic field vector Hn(r, t), is chosen. Let a unit vector in the direction of 
the chosen coordinate axis by identified as a vector e,„. A vector product, 
10 e^i^emxenCr), of the unit vector and the direction e„(r) is perpendicular to the 
direction of the vector e n (r), and hence it is contained in the plane of projection and 
is perpendicular to the coordinate axis e™. Let the magnitude of a vector obtained 
by projecting the projective component vector H m P (r-r c , 0 C , t) in the direction of the 
vector ep i, including the direction of the vector, be represented (S4) by a value of 
15 H^iV-rc, 0 C , t). That is, 

" m , 1 P (r-r c , O c , <) = H/(r-r c , 0 C , t)-e p V (4) 
e p , , = e m x e n (r). (5) 
The next step is to calculate a vector ^ perpendicular to the directions of 
the vectors ep,, and e„(r). The direction of the vector e^ is also perpendicular to 
20 the direction of the vector e„(r), and hence it is contained in the plane of projection 
and is perpendicular to the direction of the vector e^i as well. Letting the 
projection of the projective component vector H m p (r-r c , 0 c , t) in this direction be 
represented by Hm/fr-rc ©c, 0, values H,„/(r-r c , 8 C , t)ep,i and Hm/Qr-i* 0 C , fje^ 
are two independent vectors into which the projective component vector H,n P (r-r c , 
25 0 C , t) is separated. Here, 

H^ 2 v (r - r c , 0 e , ()= H ro P (r- r c , 0 C , /)• e p> 2 . (6) 



e p, 2 - e p. 1 



xe„(r> (7) 



Then, by setting the position vector r c and angle-of-orientation 0 C of the 
magnetic field source so that a theoretically calculated magnetic field, LU(r-r c , 0 C , t) 
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generated by the magnetic field source of the position vector ^substantially 
matches with a projective component He P (r-r c , 0 C , t) on the same plane as that of 
the measured magnetic field, it is possible to detect the position and orientation of 
the magnetic field source. 
5 The above description has been given of a case where the number of noise 

magnetic fields is virtually one, but when the number of noise magnetic fields is 
two, letting the direction of two noise magnetic fields H^r, t), where j = 1, 2, be 
represented by vectors eni(r) and e^O"), use is made of the projection of the 
measured magnetic field on a direction vector eN(r)=eni(r)x e^r) in the coordinate 
10 system of Fig. 3 as depicted in Fig. 5. That is, 

H m P (r-r c , 0 C9 0=H m (r-r c , t)- (H ra (r- r c , />e N (r)> N (r) (8) 

is calculated (SI 1, S12, S13). In the above 

e N (r)=e nl (r)x en2 (r) (9) 
where the symbol "x" indicates a vector product In this instance, an independent 
15 component in the projective component is one, that is, only the magnitude of the 
vector. 

Now, let the coordinates of the magnetic field source be represented by a 
vector r c (x, y, z) and its angle of orientation by a vector 0 c (0x, 0 y , 0 Z ). In the 
following description, 0^ 0 y and 0 Z will be referred to as an angle of rotation, a tilt 
20 angle and an azimuth, respectively. When the signal magnetic field is axially 
symmetric, the axis of symmetry is regarded as the x- axis, and the angle of 
orientation is set to a vector 0 c (0 y , 0 Z ). 

Further, according to the present invention, it is possible to calculate the 
magnitude H m (r-r c , 0 C ) of each component of the original signal magnetic field 
25 vector H m (r-r c , 0 C , t) by synchronous detection of the measured magnetic field 
vector Hm^-rc, 0 C , t) through the use of a proper one component in Eq. (1) or (6) 
that is a projective component. 

H s (r-r c , 0 c ) = <H m (^r c , 0 C , t)H J(r-r c9 0 c , ,)>,• (10) 

In the above, H m p (r-r c , 0 C , t) is a proper component of the projective 
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component given in Eq. (1) or (6). 

Fig. 4 shows the case where the number of noise magnetic field sources is 
essentially one. The processing of Fig. 4 is repeated at each measuring point. In 
the case of essentially two noise magnetic fields, too, the processing shown in Fig. 
5 5 is repeated at each measuring point as is the case with Fig. 4. 

The magnitude H m (r-r c , 0 C ) of the original signal magnetic field H,n(r-r c , 9 C , 
t) can be obtained by calculating Eq. (10) after calculating the projective 
component at each measuring point. Concretely, the processing depicted in Fig. 
6 or 7 is carried out [(S 11, S12, S14, S15, S16) or (Sll, S12, S14a, S16)]. 
10 How to calculate the position and angle of orientation: 

(1) Definitions of the coordinate system and the angle of orientation: 

A definition will be given first of the coordinate system necessary for the 
description of the invention. 

A coordinate system is set which is fixed to the earth with the z- axis in a 
15 vertical direction (upward), which system will hereinafter be referred to as a 
measuring coordinate system vector C M . The x- and y- axes are properly set so 
that they form a right-hand system. For example, they are set in parallel to the 
direction in which the side of a measuring frame is projected on a horizontal plane. 

This has for its object to calculate the coordinate vector r e (x, y, z) and 
20 angle-of orientation vector Q c (Qx, 0 y , 0 Z ) of the magnetic field source in the 
coordinate system. 

On the other hand, as the coordinate system vector C c of the magnetic field 
source, a coordinate system with the axial direction of the magnetic field source as 
the x axis is set so that the y- and z- axes are horizontal and vertical (upward), 
25 respectively, when the magnetic field source is placed horizontally. 

The angle-of-orientation vector 0 C of the magnetic field source is defined as 
the angle of rotation between the measuring coordinate system vector Cm and the 
coordinate system vector C c as described below. In the first place, a coordinate 
system Ceo parallel to the measuring coordinate system vector C M is turned by an 
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azimuth angle 0 Z about the z- axis (in either one of the coordinate system). This 
coordinate system will be referred to as a coordinate system vector C c i. Next, the 
coordinate system vector C cX is turned by a tilt angle 9 y about the y- axis of the 
coordinate system vector C c i itself. This coordinate system will be referred to as 
5 a coordinate system vector Further, the coordinate system vector C& is 

turned by an angle of rotation 0z about the x- axis of the coordinate system vector 
C C 2 itself. The angle-of-orientation vector 0 C is determined so that the resulting 
coordinate system becomes the coordinate system vector C c . 

(2) Explanation of independent measurands and unknowns: 

10 When the number of noise magnetic fields is essentially one, the 

measurement of the magnetic field at one place will provide two independent 
measurands. When the number of noise magnetic fields is essentially two, the 
measurement of the magnetic field at one place will provide one independent 
measurand. Further, in case of obtaining components of the original signal 

15 magnetic field by synchronous detection, the measurement of the magnetic field at 
one place will provide three independent measurands. On the other hand, three 
coordinate components of the coordinate vector r c (x, y, z) of the magnetic field 
source are unknowns. The azimuth angle 0 Z is not obtainable without a different 
criterion such as the direction of the earth magnetism, and hence it is also an 

20 unknown unless the digging head is provided with a direction sensor. When the 
digging head is made of a magnetic material, or when a buried steel pipe or similar 
magnetic body lies near the digging head, an accurate direction cannot be obtained 
even if the direction sensor is provided; hence, the azimuth angle 0 2 is an unknown 
in many cases. The tilt angle 0 y can easily be detected by a tilt angle sensor that 

25 detects the vertical direction, and hence it is known in many cases. The same is 
true of the angle of rotation 0 X . In particular, when the signal magnetic field is 
axially symmetric, if the axis of symmetry is regarded as the x- axis, the angle of 
rotation 0 X becomes meaningless and hence can be ignored. 

At any rate, measured magnetic fields need only to be obtained at different 

- i' 
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positions so that independent measurands larger in number than unknowns. 
Arrangement of the measuring system: 

For example, as shown in Fig. 1, magnetic fields are measured using a 
required number of three-axis magnetic sensors disposed on the ground so that 
5 their relative positions are known. Since the direction of the noise magnetic field 
changes for each position of measurement, it is necessary to perform for each 
position the step (S21) for detecting the position of the noise magnetic field and the 
step (S22) for calculating the projective component of the measured magnetic field 
as shown in Fig. 8 that is a flowchart of the procedure for position determination 
10 by the present invention. In case of calculating the component of the original 
signal magnetic field by synchronous detection, the processing therefor (in a 
specified step) needs to be carried out for each position. 

Flow of measurement processing: 

A detailed description will be given of the detection (S21) of the direction 
15 of the noise magnetic field, the calculation (S23) of the position of the signal 
magnetic field source through the use of the projective component and the 
calculation of the position f of the signal magnetic field through the use of the 
signal magnetic field component. 
[EMBODIMENTS] 
20 Embodiment in the case of virtually one noise magnetic field: 

In case of using the projective component, when the number of unknowns 
is Nu(>l) and the number of noise magnetic fields is virtually one, the projective 
component vector H m p (r-r c , G c , t) expressed by Eq. (1) is calculated from magnetic 
vectors H m (r-r c , 0 C , t) measured at Ny/2 or more different positions, and the 
25 position vector r c and the angle-of-orientation vector 0 C are determined (S31) as 
depicted in Fig. 9 so that a projective component vector at each magnetic field 
sensing point is essentially equal to the projective component vector He P (r-r c , 0 C , t) 
of a theoretically calculated magnetic field of the projective component at each 
position of measurement, by which it is possible to detect (S32) the position and 
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orientation of the magnetic field source. 

In case of using the magnitude of the signal magnetic field synchronously 
detected by the projective component, the magnetic field is measured at Nu/3 or 
more different positions. When the number Nu(>l) of unknowns is even, the 
5 number of unknowns and the number of independent measurands can be made 
equal to each other; if setting magnetic field vector H^r-rc, 0 C , t) is measured at 
positions N m =Nu/2, a Nu number of such equations as given below need only to be 
solved. 

K/fe"*. O e ) = 0,k = l,...,N m ,g = l,2. (11) 

10 where q = 1, 2 and represents two directions parallel to the plane of projection but 
parallel to each other (S33). Accordingly, H^V-To 0 to t) and He, q P (r-r c , 9 C , t), 
where q = 1, 2, are the magnitudes of q-direction components of the measured 
magnetic field vector H m (r-r c , 0 C , t) and the theoretical calculated magnetic field 
vector He(r-r c , 0 C , t), respectively. Let it be assumed that the vector 9 C represents 

1 5 any one of 9*, 0 y , 9 a (9 y , 9 Z ), (0 a 9 X ), (Ox, 9 y ), (9*, 0 y , 9 Z ) and (J) and that the vector r c 
represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) and <J>, where <|> represents 
an empty set 

For example, when unknowns are the position vector r c (x, y, z) and 
azimuth angle 9 Z of the magnetic field source, the magnetic field is to be measured 
20 at two different positions and four equations such as given below are solved, by 
which it is possible to obtain (S34) the position vector r c (x, y, z) and azimuth angle 
9 Z of the magnetic field source. 

K /fe - r e , 0 Z , 0), - H % q P (r k - r c , 0 Z ) = 0, * = 1, 2; q = 1, 2. (12) 

In the above, < . > t represents a time average. 
25 As depicted in Fig. 10, when the magnetic field is measured at N m (>Ny/2) 

different places more than Ny/2 in excess of the number Nu of unknowns, since a 
larger number of independent measurands than the number of unknowns can be 
obtained (S41, S42, S43), the position r c (e, y, z) and the angle of orientation 0 Z are 
calculated which provide 
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™n{||X ,|K /fe - r c , *. /ft - q P (r k - r c , * c )|. (13) 

where < . > t represents the time average and min { .} means that vectors r c and 0 C 

are changed to obtain vectors r c and 0 C that minimize the contents in { . }. Further, 
a symbol w M indicates weighting. The Eq. (13) can be replaced with the 
5 following equations (S44). 

£in{!|>, \j(K;(r k -r c ,0 c , t f\ - //„ q P (r t - r c , * c |. (14) 

^{||-,jK<, P ( r k - r c. ^/ft-^„ P (r k -r c ,^}. (15) 

min{|| Wt ,^(K q P (r k -r c , 0„ tf) - ,'(r fc - r c , ^ 2 ||. (16) 

Assume that vectors r c and 9 C in Eqs. (13), (14), (15) and (16) have the same 
1 0 meaning as in the case of Eq. ( 1 1 ). 

In the above description it does not matter whether the measured magnetic 
field vector H m (r-r c , 0 C , t) is a signal having passed through a band pass filter that 
permits the passage therethrough of only components close to the frequency of the 
signal magnetic field or a wide-band signal that is inhibited from the passage 
15 through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of deteimining the position of the 
magnetic field with high reliability. 

(Embodiment in the case of essentially two noise magnetic fields) 
When the number of unknowns is Nu (>1) and the number of noise 
20 magnetic fields is virtually two, the projective component vector H m p (r-r c , 0 C , t) 
expressed by Eq. (4) is calculated from magnetic vectors H,n(r-r c , 0 C , t) measured 
at Nu/2 or more different positions, and the position vector r c and the 
angle^of-orientation vector 0 C are determined (S51, S52) as depicted in Fig. 11 so 
that the above-mentioned projective component vector obtained at each position of 
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measurement essentially matches with the projective component vector He P (r-r c , 9 C , 
t) of a theoretically calculated magnetic field at each position of measurement, by 
which it is possible to detect (S53) the position and orientation of the magnetic 
field source. 

In this case, the number Nu (>1) of unknowns and the number of 
independent measurands can be made to be equal to each other without fail; the 
magnetic field vector H m (r-r c , 9 C , t) is measured at different positions of the same 
number as that N v number of unknowns, and Nu number of such equations given 
below need only to be solved. 

(// m P (r k -r c , 0 c ,')),-"e P (r k -r c , 0 C )=O, k = \, ... ,N„. (17) 

In this case, symbols H m p (r-r c , G c , t) and He P (r-r c , 9 C , t) are the respective 
magnitudes of the projective component vector H m p (r-r c , 9 C , t) and the vector 
H/Cr-rc, 0 C , t) of the measured magnetic field vector HJ{r-r c , 9 C , t) and the 
theoretical calculated magnetic field vector He(r-r c , 9 C , t). Let it be assumed that 
the vector 9 C represents any one of 9 to 9 y , 9 a (9 y , 9 Z ), (9 & 9*), (9*, 9 y ), (9*, 9 y , 9 Z ) 
and (j) and that the vector r c represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) 
and (J), where <|) represents an empty set. 

For example, when unknowns are the position vector r c (x, y, z), azimuth 
angle 9 Z and tilt angle 9 y of the magnetic field source, the magnetic field is to be 
measured at five different positions and four such equations given below are 
solved, by which it is possible to obtain the position vector r c (x, y, z), azimuth 
angle 9 Z and tilt angle 9 y of the magnetic field source. 

(tf ra P (r k - r c , <9 c ,/)),-// e P (r k -r c , 0 C )=O, k = \, ... ,5. (18) 

In the above, < . > t represents a time average. And the vector 9 C = 9 c (9 y , 9 Z ). 

As depicted in Fig. 12, when the magnetic field is measured at N m (>Nu) 
different places more than Nu in excess of the number Nu of unknowns, since a 
larger number of independent measurands than the number of unknowns can be 
obtained (S61, S62), the position r c (e, y, z) and the angle of orientation 9 Z are 
calculated which provide 
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min< 




(19) 



where < . > t represents the time average and min { .} means that vectors r c and 9 C 

r c»°z 

are changed to obtain r c and 0 C that minimize the contents in { . }. Further, a 
symbol w^q indicates weighting. The Eq. (19) can be replaced with the following 
equations (S63). 



Assume that r c and 0 C in Eqs. (19), (20), (21) and (22) have the same meaning as in 
the case of Eq. (16). 

In the above description, it does not matter whether the measured magnetic 
field vector H m (r-r c , 0 C , t) is a signal having passed through a band pass filter that 
permits the passage therethrough of only components close to the frequency of the 
signal magnetic field, or a wide-band signal that is inhibited from the passage 
through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of determining the position of the 
magnetic field with high reliability. 

Further, in case of using the signal magnetic field component obtained by 
synchronous detection, the direction of the noise magnetic field is determined and 
is used to obtain the projective magnetic field; and the subsequent processing is 
common to the cases of one and two noise magnetic fields. The projective 
component vectors H^r-rc, 9 C , t) at respective places are calculated from 
magnetic field vectors H m (r-r c , 0 C , t) measured at Ny/3 or more different places, 
and a proper one of the calculated vectors is used as a reference signal to perform 




(22) 



(20) 



(21) 
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synchronous detection of the measured magnetic field vectors H m (r-r c , 0 C , t), 
thereby obtaining the magnitude Hs(r-r c , 0 C ) of the original signal magnetic field 
component. By determining the position vector r c and the angle-of-orientation 
vector 0 C of the magnetic field source so that the magnitude of the original signal 
5 magnetic field component and the magnitude Hc(r-r c , 9 C ) of the theoretical signal 
component are equal to each other, it is possible to detect the position and 
orientation of the magnetic field source. 

When the number Nu (>1) of unknowns is a multiple of 3, the number of 
unknowns and the number of independent measurands can be made to be equal to 
10 each other; if the magnetic field vector H m (r-r c , 0 C , t) is measured at positions 
Nm^Nm/Nu/3, Ny number of such equations as given below need only to be solved. 

H s (l-r c , 0>H c (r k -r c , 0 c )= 0 ,* = 1, (23) 

Fig. 13 shows the flow of processing (S61, S62, S65). Here, H s (r-~r c , 0 C ) and 
Hc(r-r c , 0 C , t) are the signal magnetic field calculated from the measured magnetic 

15 field vectors H m (r-r c , 0 C , t) and the magnitude of the theoretical calculated 
magnetic field vector, respectively. The magnitude Hs(r--r c , 0 C ) of the signal 
magnetic field vector is an averaged quantity as already explained. Let it be 
assumed that the vector 0 C represents any one of 0^ 0 y , 0 to (0 y , 0 Z ), (02, 0 X ), (0x, 0 y ), 
(0x, 0 y , 0 Z ) and <|> and that the vector r c represents any one of x, y, z, (y, z), (z, x), (x, 

20 y), (x, y, z) and (J>, where (J) represents an empty set. 

For example, when an unknown is the position vector r c (x, y, z) of the 
magnetic field source, the magnetic field is to be measured at one place and three 
such equations as given below are solved, by which it is possible to obtain the 
position vector r c (x, y, z) and azimuth angle 0 Z of the magnetic field source. 

25 H s (r-r c )-H c (r-r c )=0. (24) 

The vector r is the vector of the position of measurement. 

When the magnetic field is measured at N m (>Ny/2) different places more 
than Ny in excess of the number Nu of unknowns, since a larger number of 
independent measurands than the number of unknowns can be obtained, the 
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position r c (e, y, z) and the angle of orientation 0 Z are calculated which provide 
min{i>jH,0- k " * c > *.)" H c (r k - r c , *jl (25) 

where min { . } means that vectors r c and 0 C are changed to obtain vectors r c and 0 C 

that minimize the contents in { . }. Further, a symbol w^q indicates weighting. 
5 The Eq. (25) can be replaced with the following equations. 

^in{Z^(|Hs(r k - r c , *.| - |H.(r k - r c , 0 C |) 2 J. (26) 

min||j^|H i fe - r c , *> « c (r k - r c , (27) 

Fig. 14 shows the flow of processing (S61, S64, S66). 

A description will be given below of how to determine the direction of the 
10 noise magnetic field in the two embodiments described above. 

(How to determine the direction of the noise magnetic field in the case of 
virtually one noise magnetic field) 

(First method) 

The first method for determining the direction of the noise magnetic field is 
15 a method where, in case of no signal magnetic field, a noise magnetic field is 
measured through the use of the same measuring system as that for measuring the 
signal magnetic field. This situation is such as shown in Fig. 15 in which the 
detection of the digging position is disturbed by a noise magnetic field source lying 
in the vicinity of the digging route. Another case is that where the signal 
20 magnetic field source is equipped with a function of receiving a command sent, for 
example, from the ground by some means and responsive to the command to stop 
the generation of the signal magnetic field. 

In this instance, as depicted in Fig. 16, letting the measured magnetic field 
be represented by a vector H m (r-r c , 6 C , t), since this is essentially a noise magnetic 
25 field vector H„(r, t) as shown at a step (S71), an average value of its absolute values 
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can be used (S72) to calculate a direction vector, en(r)=en, x (r), en, y (r), e^r)), of the 
noise magnetic field as follows: 

e n . a 00 = -7— = ^ , a=x,y, z. (28) 

Alternatively, a root-mean-square value of the above absolute values can be used to 
determine the direction of the noise magnetic field by 



, , a = x,y,z. (29) 

\a = x,y,z 

In this case, a symbol H m ^(r-r c , 9 C , t) is an a component of the measured magnetic 
field (the noise magnetic field), and a is any one of x, y and z. 

In the above description it does not matter whether the measured magnetic 
10 field vector H m (r-r c , 0 C , t) is a signal having passed through a band pass filter that 
permits the passage therethrough of only components close to the frequency of the 
signal magnetic field or a wide-band signal that is inhibited from the passage 
through the band pass filter, but the use of the signal having passed through the 
band pass filter increases the possibility of determining the position of the 
1 5 magnetic field with high reliability. 
(Second method) 
Step 1 

As shown in Fig. 17, at first , the frequency spectrum H m (co) of the 
measured magnetic field vector H m (r-r c , 0 C , t) is calculated (S81, S82, S83) by the 
20 following equation. 

H m (a>)=F(H m (r k -r c , /)) (30) 
In this case, a symbol F( . ) represents a Fourier transform; the three components x, 
y and z of the measured magnetic field vector H m (r-r c , 0 C , t) are each Fourier 
transformed. In practice, the above-mentioned frequency spectrum can be 
25 calculated by FFT (fast Fourier transform) or the like of sampled values of the 
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measured magnetic field vector H m (r-r c , 0 C , t). 
Step 2 

Next, an angular frequency ©i, where i = 1, 2, N s , of a large-amplitude 
component, such as a line spectrum, is selected (S84) from the absolute value 

5 |H m ((o)| of the frequency spectrum. For the component of each angular frequency 
©i, where i = 1 , 2, . . . , Nns , a candidate unit vector e„(r, ©0, where i = 1, 2, . . . , N^, of 
the direction of the noise magnetic field, is calculated (S85) by the method (1) or 
(2) described below. 

(1) The absolute values of the Fourier-transformed x, y, and z components 

10 of the angular frequency concerned are used to calculate the candidate unit vector 
e„(r, ojj), where i == 1, 2, . . . , N„s, of the direction of the noise magnetic field, by the 
following procedures: 

( x (H m Jfi4 |H m ,Xo>$ |H ' 
e a (r,a> i )= ■ , „ =g > ' = 1, ... , N m . (51) 



\a=x,y,z 



where H^Oi), where a = x, y, z and i = 1, 2, . . . , N„s, is a ©i component by the 
1 5 Fourier transform of the a component of the measured magnetic field. 

(2) A narrow-band filter is formed whose pass band uses, as the center 
frequency, the angular frequency cOj, where i = 1, 2, Nns , and the same method 
as by Eq. (21) or (22) is used to calculate the candidate unit vector e„(r, ©;), where i 
= 1, 2, . . . , N„s . That is, the candidate unit vector e„(r, a>i)= (e^O"), e„, y (r), e,^(r)), 
20 where i = 1, 2, . . . , N m ,is calculated by 

JZ K .( r - r «- <>.•<».' 'Hi 

f a=jr, y, z 

or by 



(r, 6>j)= i , , a = .V, 2; ' = 1, ■• , 

J Z <(^. a (r k -r c , «.«k./)y) 

I a~x, y, z 



(33) 
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Step 3 

The candidate unit vector e„(r, ©i)= (e^r), e„j<r), e„^(r)), where i = 1, 2, . . . , 
Nns , is considered as the direction vector e„(r) of the noise magnetic field, and for 
each angular frequency ©j, where i = 1, 2, . . . , N„s , the same method as by Eq. (1) 

5 is used to calculate a projective component vector H m p (r-rc, 0 C , C0j, t). 
H m P (r-r c , 0 C , a; i5 0=H m (r-r c , 0 C , t) 

-(H m (r-r c , 0 C , f).e.(r. ^)> n (r, «J i = h ... >„. (34) 
The reason for which the angular frequency ©i is contained as a variable of the 
projective component is to explicitly point out that the projective component is 

10 dependent on the angular frequency ©j, where i = 1, 2, N„s . A roper time 
interval T^, which consists of durations T^, where k = 1, 2, N^, each 
having a short time length At, is chosen, and the variation of the projective 
component vector H m p (r-r c , 0 C , ©;, t) for each duration T^ where k = 1, 2, . . . , Ntest, 
is evaluated (S86). Assume, here, that each duration T^ where k = 1, 2, . . . , 

15 does not overlap other durations. Concretely, a variance of statistics v^^©;) 
of the N^, where k = 1, . . . , N test , which are calculated by any one of the methods 
described below, is calculated. 

(1) One or both of the means of absolute values of two orthogonal components by 

v e v al , M> (K, /fr " r c . o>, k , q = 1, 2; k = 1, ... , A^; / = 1, ... , N n , (35) 

20 where ( ) Tt k represents the mean value in the duration Tgc and VevaUc is a statistic 
calculated for the duration T^. 

(2) Mean of absolute values 

Veva,, k («i)=(|H m P (r-r c , 3, co„ ty^, k = \ ... , N^; i= 1, ... , N m . (36) 

(3) One or both of the means of squares of two orthogonal components by 
25 k (*0=((/^ q P (r- r c , 0 C , a>„ ')j) T L 2; *- 1, ... , i= 1, ... , N m . 

(37) 
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(4) One or both of square roots of the means of squares of two orthogonal 
components by 



V . 
cval, 



k «= J{K:(r- r *> 0l\ k - 9 = 1. 2; k = 1, ... , A- i = 1, ... , N m 



(38) 

5 For the statistics Ve^, where k = 1, calculated by these equations, the 

following equation 



Jmean k (fv^ k (a>,y mean^v^ k (<yj))) 2> ) 
var(o> i )= J * ^ -r-^r 1, ... . (39) 

is calculated to obtain (S86) a value of m i/nin that is the angular frequency a>; which 
minimizes var(a>i). In the above, mean k (.) indicates averaging for the suffix k, 
10 that is, 

So 

mean k (.) = ^- . (40) 

The magnetic field of the angular frequency derives from the noise 
magnetic field, and the direction of the noise magnetic field becomes a vector e„(r, 

®i^nin)- 

15 Incidentally, the angular frequency a),^ which minimizes var(G>i), needs 

only to be measured at one place and need not be obtained at every place where to 
measure the magnetic field. 

With this method, it is also possible to calculate a fluctuation in the 
direction of a vector H m p (r-r c , 9c, ©i, t) as well as an amplitude fluctuation given by 
20 Eq. (39) and to select the angular frequency co^ at which the direction fluctuation 
becomes minimum or smaller than a predetermined value. 

Incidentally, in the above description, the measured magnetic field vector 
H m (r-r c , 9 C , t) in Step 1 is a wide-band signal, and in Step 3 it does not matter 
whether the measured magnetic field vector H m (r-r c , 0 C , t) is a signal having passed 
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through a band pass filter that permits the passage therethrough of only 
components close to the frequency of the signal magnetic field or a wide-band 
signal that is inhibited from the passage through the band pass filter, but the use of 
the signal having passed through the band pass filter increases the possibility of 
5 determining the position of the magnetic field with high reliability. 

Fig. 18 shows how to select frequencies fi(=coi/27c), f 2 (=<02^2n% 
f n (=(G>ii/27r, where n=N ns ), at which the frequency spectrum becomes maximum; 
Fig. 19 shows the flow of processing including step S91 for obtaining the 
candidate vector; and Figs. 20 and 21 show the flow of processing including steps 
10 S101 and S102 or steps Sill and SI 12 for evaluating the candidate vector and for 
detecting the direction of the noise magnetic field. 

The frequency spectrum H^co) need not always be used. That is, the 
signal magnetic field is periodically turned OFF/ON following a predetermined 
procedure; the period Tp^ is divided into equally-spaced durations; the candidate 
1 5 unit vector e„(r, ti), where i = 1 , . . . , N,*, is used in place of the candidate unit vector 
e n (r, coO, where i = 1, . . . , which is calculated by Eqs. (25), (26) and (27); and 
thereafter, the duration that minimized the variance by Eq. (33) is calculated by the 
processing described above. By this, the vector en(r, ti) in that duration can be 
adopted as the direction of he noise magnetic field. 
20 (Third method) 

In the second method, the large-amplitude angular frequency where i = 
1, 2, . . . , N s , is selected from the absolute value |H m (co)| of the frequency spectrum 
H m (co), but the vector e^r, oo^) can be obtained as the direction of the noise 
magnetic field in exactly the same manner as in the case of the second method, by 
25 selecting a proper frequency band neighboring the frequency of the signal 
magnetic field, setting properly-spaced test frequencies free from the frequency of 
the signal magnetic field in the frequency band and regarding the test frequencies 
as the angular frequency (0\ in the second method. 

As is the case with the second method, the angular frequency at 
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which var(cDi) becomes minimum needs only to be obtained. 

Fig. 22 shows how to select the frequencies fi(=coi/27c), f 2 (=G)2/27c), 
fnCK^n^rc, where n=N lls ) at which the frequency spectrum becomes maximum. 
The flow of the subsequent processing is the same as depicted in Figs. 17, 19 and 
5 20. 

(Fourth method) 

The signal magnetic field is periodically stopped under control of a 
predetermined procedure. For example, the signal magnetic field is periodically 
stopped by a predetermined time interval. Since the intensity of the magnetic 

10 field being measured decreases while the signal magnetic field is stopped, the OFF 
period of the signal magnetic field is identified by regarding the 
intensity-decreasing period essentially as the predetermined OFF period, and the 
direction of the magnetic field measured during the OFF period is used as the 
direction of the noise magnetic field. The direction of the noise magnetic field 

15 can be obtained using the same method as the first one. Fig. 23 shows how the 
amplitude of the measured magnetic field in this method varies with time. 
(Fifth method) 

The signal magnetic field is periodically stopped following a 
predetermined procedure. This is carried out as described below under (1) and 

20 (2). 

(1 ) To stop the signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the period 
Tpcnod, for instance. 

= <^<«,+W ( 

25 and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 

the sequence is -1, the signal magnetic field is turned ON. In the above, 

'w-' k = W *=1>2, 3,... (42) 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise 
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basis: 

For example, when the value is in a random sequence like an 

M-sequence consisting of unit periods of the same length Nm, the signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
5 turned OFF; and this sequence is repeated. In this case, the time average of the 

sequence is set to 0. 

Fig. 24 shows temporal variations in the amplitude of the measured 
magnetic field in the case of the signal magnetic field being stopped by the method 
(1). Fig. 25 shows temporal variations in the amplitude of the measured magnetic 
10 field in the case of the signal magnetic field being stopped by the method (2). 

Next, as depicted in Fig. 26, the correlation function between the sequence 
s(t) and the norm of the measured magnetic field or the absolute value of its 
particular component is calculated (S121, S122, S123). As the correlation 
function, it is possible to use any one of those calculated by 

15 1<t) = r^lH^r-r., 0 C , ifc-vy*. (43) 



KrH W ^V(">- r <> * rH (44) 

ft**'" K air - r c , 4, 4* - r)dt, a = x,y,z. (45) 
^ a(r _ re ,^ t )f <t - T y t , a = x,y,z. (46) 

In this case, the period for which to detect the correlation is set to an integral 
20 multiple N T Tperiod of the period 

Tperiod- 

The OFF state of the signal magnetic field can be detected (SI 25) from the 
time T=t sync at which any one of the above correlation functions become maximum 
(S 1 24). That is, a sequence s(t-tsyncX which starts at the time T=t sync , is used, and 
when the sequence s(t) is "1," the signal magnetic field is regarded as being OFF; 
25 in this way, the ON/OFF operation of the signal magnetic field is determined. 
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In the thus determined signal magnetic field OFF period the direction of 
the measured magnetic field is dtetefited, and the direction is regarded as the 
direction of the noise magnetic field (S126). When the sequence s(t) is (a 
second numerical value), the signal magnetic field is turned ON, and when the 
5 sequence is "1" (a first numerical value), the signal magnetic field is turned OFF; 
this sequence is repeated in this way. 
(Sixth method) 

As is the case with the fifth method, the signal magnetic field is 
periodically stopped under control of such a predetermined procedure as described 
10 below. 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the 
period T^od, for instance. 

s(t) = l, +W ( 41 ) 

= -1, t k +t Aop <t<t k ^. 

15 and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
the sequence is - 1 , the signal magnetic field is turned ON. In the above, 

' k+1 -<* = W*=L2, 3,... < 42 > 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise 

basis: 

20 For example, when the value is in a random sequence like an 

M-sequence consisting of unit periods T^ of the same length N M , the signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
turned OFF; and this sequence is repeated. 

In this case, the sequence s(t) is chosen so that it changes for each 
25 predetermined time unit AW- And, the time average of the sequence is "0." 

Next, as depicted in Fig. 27, the correlation function between the sequence 
s(t) and the norm of the measured magnetic field or the absolute value of its 
particular component is calculated (SI 31, S132, S133) as in case of the fifth 
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method. As the correlation function, it is possible to use any one of those 
calculated by 

R(t) = f, W "> ra (r - r c , 0 C , ijpd - t}*. (43) 
Kr) = j;^^( Hm(r _ rc) I ,))%(, _ r y, (44) 

5 ^(r)= j'^^ K «(■* " <K' " r >"» « = W- (45) 

^(r)= £ +A ^ V(^-( r - r c»^ 0^/ - r)4, « = x,*z. (46) 

In this case, the period for which to detect the correlation is set to an integral 
multiple NxTpenod of the period T^od 

In this instance, there are present, in general, plural times T=t synCjk (k = 1, 2, 
10 . . . , Ngync) in which the correlation function becomes maximum and the maximum 
value exceeds a predetermined value (SI 34). Assume, for example, that t^k, (k 
= 1, 2, N^) is an arrangement of such times in order of time. When the 
correlation value between the sequence s(t) and the signal magnetic field is 
appropriate, 

1^ A*sync, k = 'sync, k ~~ 'sync, 1 > & = 2, ...» ^sync (^7) 

is virtually an integral multiple of the time unit At^ Then, the average of the 
value resulting from the subtraction of an integral multiple M^^At^ of the time 
unit Aty™, where k = 2, ^ from At^k , where k = 2, N^, is 

calculated by 

2 (A/ sync , k - Mync, k A/ unit ) 

20 <* sync =^ (48) 

iV sync 1 



In this case, 

'sync = ^sync, 1 ^sync (^9) 

provides the beginning of the sequence signal corresponding to the ON/OFF 
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operation of the signal magnetic field. 

Accordingly, the period during which the signal magnetic field is OFF can 
easily be set based on the sequence s(t-tsync)- 

By applying the same method as the first method to the magnetic field 
5 vector H m (r-r c , 0 C , t) measured in this period, the direction e„(r) of the noise 
magnetic field vector Hn(r, t) can be calculated (S 138). 

(Seventh method) 

This method will be described below with reference to Figs. 28, 29, 30 and 

31. 

10 As is the case with the fifth method, the signal magnetic field is 

periodically stopped, for example, by such a procedure as described below. 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the 
period T^od, for instance. 

s(0 = h t k <i <t k +/ t . 
15 stop (41) 

= -l> '* + '.op ^ 'inl- 
and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
the sequence is -1, the signal magnetic field is turned ON. In the above, 

'k+l ~~'k = period ^ h 2, 3, ... (42) 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise 

20 basis: 

For example, when the value is in a random sequence like an 

M-sequence consisting of unit periods T^ of the same length Nm, the signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
turned OFF; and this sequence is repeated accordingly. In this case, the time 
25 average of the sequence is "0." 

Next, the correlation function between the sequence s(t) and the measured 
magnetic field HUr--r c , 0 c , t) is calculated (S141, S142, S143). The period Tp^od 
is divided into equally spaced N^ sections of a length T^y, and either one of the 
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following calculations is conducted. 



^(0=C +rpatod K«(r-r c , 4.t}(i-i k y*.k-l....N m ,ia-xw. (50) 

R*(ft)= £ ^jQtm. a(r - r c , 4, O) 2 «0 " 'kK *=L ... »^<«v; « = (51) 

In this case, a symbol Ra(tk), where a = x, y, z and k = 1, N^, is the time 
correlation between an a component H^r -r c , 0 C , t) of the measured magnetic 
field H m (r-r c , 0 C , t) and the sequence s(t). And 

t k =t 0 +k-T^ k = l...,N div . (52) 
The measured magnetic field may also be correlated with the time that is an m- 
multiple of the period T^od- That is, 



10 4. '&,('- 0<*. * =1 » « = W- (53) 

^ (/) = £~-fa m (r-^ 0 e , iffs^Q-tjik. * = ...,N div ; a=x,y,, (54) 
In this instance, the sequence s(t) is replaced with Smp(t), where 

^pCO^^+^od) (55) 

S mp (t)=s(0, o<t<r period . 

This is followed by calculating the component H^ir-Tc, 0 C , tk, t), where k = 
15 1 , . . . , Ndiv , of the measured magnetic field H m (r-r c , 0 C , t) projected on the vector 
e n (tkMR x (t k ), Ry(tk), IW), where k = 1, . . . , 
formed by correlation fiinctions R«(tkX where a = x, y, z , corresponding to the 
respective components x, y and z of the measured magnetic field H m (r-r c , 0 C , t). 
Here, a symbol tk contained as a variable of the projective component indicates that 
20 the projective component depends on the variable tk. 

In this method, the vector e^k) can be used as the direction of the noise 
magnetic field through utilization of the time tk in which a fluctuation in the 
absolute value of the projective component vector H^r-rc, 0 C , tk, t), where k = 1, 
. . . , Ndi v , 
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var 



J((HL P (r-ii, Z, 0-(|H m p (r-r c , t k , 

var(, k )=-U 7 — — — =- ± * = 1, . ^ div . (56) 

flB„ (r-r c , O c , t k , tj[) 

becomes minimum or smaller than a predetermined value (SI 45a). In the above 
() t means the calculation of the time average. 

Further, variances of the x-component H^r-rc, 9 C , tk, t), y-component 
Hn^y^-rc, 0 C , tk, t) and z-component H^r-rc, 0 C , tk, t) of the projective component 
vector H m (r-r c , 0 C , tk, t), where k = 1 , . . . , N^, 

_ ^((" m ,/(r-r 0 , 0, Z 0-K .'('-«•-. *o Z 'J>J) 

a = x,y 9 r, * = 1, ... , Af div . 

are calculated, and the vector e n (tk) can be used as the direction of the noise 
10 magnetic field through utilization of the time tk in which the sum of the 
above-mentioned variances 

Z var„(0 (58) 

a=x. y. z 

or 

J Z (varAtf (59) 

15 becomes minimum or smaller than a predetermined value (S145b). In this case, 
there is the possibility that the correlation functions R x (tk), R y (tk) and R z (tk) at a 
certain time tk have lost their original signs. Hence, it is necessary to evaluate the 
fluctuation of the projective component at each time tk for four combinations 

[RxOk), R y (U Rz(tk)L [R x (tkX R y (tk), -R z (tk)L [Rx(U -RyO*), R 2 (tk)] and [R x (tk), 
20 -R y (U -R z (tk)]. 

Further, the period during which the signal magnetic field is OFF can 
easily be set (S145c, S145d) based on the sequence s(t-tk). By applying the same 
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method as the first one to the magnetic field vector H n (r-r c , 0 C , t) measured in this 
period, the direction e^r ) of the noise magnetic field vector H n (r, t) can be 
obtained (SI 46). In the above, use is made of the time k of the equally spaced 
Ndiv sections of a the time length divided from the period Tp^, but it is also 
5 possible to use the time when the correlation function given by Eq. (45) or (46) 
^ becomes maximum or when the correlation function becomes maximum and 

exceeds a predetermined value. 

(Method for determining the direction of the noise magnetic field when the 
number of noise magnetic fields is virtually two) 
10 Even in a case of two noise magnetic fields, when a first one of the two 

noise magnetic fields has far higher an intensity than the second noise magnetic 
field at a first frequency and the second noise magnetic field has far higher an 
intensity than the first noise magnetic field at the second frequency, the directions 
of the first and second noise magnetic field can easily be calculated through 
15 utilization of these frequency components in the measured magnetic field. 
^ If the first or second frequency is close to the frequency of the signal 

magnetic field, the directions of the noise magnetic fields can be determined from 
the measured magnetic field having passed through a band pass filter that permits 
the passage therethrough of only frequencies near those of the signal magnetic field, 
20 by the same method as the fourth or fifth method for use in the case of virtually one 
noise magnetic field 

When either of the first and second frequencies does not equal to the 
frequency of the signal magnetic field, the directions of the respective noise 
magnetic fields need only to be calculated by the same method as the first one for 
25 use in the case of virtually one noise magnetic field 
(Other embodiments) 

The present invention is also effective when the signal magnetic field 
generated by the magnetic field source is virtually axially symmetric, and the 
invention permits position determination with a smaller number of magnetic 
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sensors or by magnetic field sensing at a smaller number of positions than in the 
case of a magnetic field of low symmetry. 

Further, when only one noise magnetic field affects the measurement of the 
digging position and the tilt angle of the magnetic field source is known which is 
5 an inclination of the axis of symmetry corresponding to the axial direction of the 
signal magnetic field set in die magnetic field source with respect to the vertical 
direction, the projective component of the magnetic field, measured at each of two 
or more different positions, on a plane perpendicular to the direction of the noise 
magnetic field sensed at each of the magnetic field sensing positions is calculated; 

10 the position of the magnetic field source and its azimuth angle that is the direction 
of the axis of symmetry in a horizontal plane can be calculated from the above-said 
projective component. 

Further, according to this invention method, when virtually one noise 
magnetic field affects the measurement of the digging position, the projective 

15 component of the magnetic field, measured at each of three or more different 
positions, on a plane perpendicular to the direction of the noise magnetic field 
sensed at each of the magnetic field sensing positions is calculated; the position of 
the magnetic field source, its tilt angle that is an inclination of the axis of symmetry 
corresponding to the axial direction of the signal magnetic field set in the magnetic 

20 field source with respect to the vertical direction, and the azimuth angle of the 
magnetic field source that is he direction of the axis of symmetry in the horizontal 
plane can be calculated from the above-said projective component. 

Further, when virtually two noise magnetic fields alone affect the 
measurement of the digging position and the tilt angle of the magnetic field source 

25 is known which is an inclination of the axis of symmetry corresponding to the axial 
direction of the signal magnetic field set in the magnetic field source with respect 
to the vertical direction, the projective component of the magnetic field, measured 
at each of four or more different positions, on a straight line perpendicular to both 
of the direction of a first one of the two noise magnetic fields sensed at each 
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magnetic field sensing position and the direction of the remaining second noise 
magnetic field sensed at the same position is calculated; the position of the 
magnetic field source and its azimuth angle that is the direction of the axis of 
symmetry in the horizontal plane can be calculated from the above-said projective 
5 component. 

|p Further, when virtually two noise magnetic fields alone affect the 

measurement of the digging position, the projective component of the magnetic 
field, measured at each of five or more different positions, on a straight line 
perpendicular to both of the direction of a first one of the two noise magnetic fields 
10 sensed at each magnetic field sensing position and the direction of the remaining 
second noise magnetic field sensed at the same position is calculated; the position 
of the magnetic field source, its tilt angle that is an inclination of the axis of 
symmetry corresponding to the axial direction of the signal magnetic field set in 
the magnetic field source with respect to the vertical direction, and the azimuth 
15 angle that is the direction of the axis of symmetry in the horizontal plane can be 
^ calculated from the above-said projective component. 

Moreover, when virtually two noise magnetic fields alone affect the 
measurement of the digging position, the frequency component of a first one of the 
two noise magnetic fields, in the vicinity of which the remaining second noise 
20 magnetic field and the signal magnetic field have substantially no frequency 
components, is measured to thereby permit detection of the direction of the first 
noise magnetic field in terms of vector; and the frequency component of the 
second noise magnetic fields, in the vicinity of which the first noise magnetic field 
and the signal magnetic field have substantially no frequency components, is 
25 measured to thereby permit detection of the direction of the second noise magnetic 
field in terms of vector. 

In the present invention, it is effective to use, as a magnetic sensor, a 
three-axis magnetic sensor that senses three magnetic fields orthogonal to one 
another at substantially the same position. 
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The magnetic sensor for use in the present invention may be any kinds of 
sensors as long as they are capable to three magnetic fields orthogonal to one 
another at substantially the same position, but the three-axis magnetic sensor is 
suitable which senses three magnetic fields orthogonal to one another at 
5 substantially the same position. Alternatively, it is possible that one magnetic 
sensor capable of sensing a magnetic field in only one direction is turned at the 
same position toward three orthogonal directions one after another to sense the 
three orthogonal magnetic fields. 

In carrying out the present invention, it is possible to employ such a frame 

10 12 as depicted in Fig. 32 which has magnetic sensor fixing means 11 mounted 
thereon to fix three-axis magnetic sensors and a tilt angle gauge 13 for detecting 
the inclination of the frame with respect to a vertical direction. The position of 
each magnetic sensor fixing means on the frame is known, and the magnetic sensor 
fixing means possesses a function of fixing the magnetic sensor in a predetermined 

15 orientation to the frame. The magnetic sensor fixing means 11 is provided with, 
for example, three faces orthogonal to one another, and has a mechanism that fixes 
the magnetic sensor at a predetermined angle when a predetermined face of the 
sensor case is pressed against any one of the three faces of the sensor fixing means. 
One magnetic sensor or one three-axis magnetic sensor is fixed to these magnetic 

20 sensor fixing means one after another to sense the magnetic fields. 

In another alternative, a plurality of magnetic sensors may be fixed in 
predetermined orientations to the frame 12 at a plurality of positions to 
simultaneously sense magnetic fields at the plurality of positions. 

As described above, the present invention employs a frame provided with a 

25 plurality of magnetic sensor fixing means each capable of removably or fixedly 
mounting a three-axis magnetic sensor and a tilt angle sensor capable of detecting 
the tilt angle of an orthogonal coordinate system to the vertical direction, the 
magnetic fixing means being mounted on the frame so that their positions and 
orientations are known; the magnetic sensor is removably or fixedly mounted on a 
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required number of magnetic sensor fixing means to sense magnetic fields, and the 
tilt angle of the frame during magnetic field sensing and the orientation of the 
magnetic sensor at each magnetic sensor mounting position with respect to the 
frame are used to calculate from the magnetic field sensed at each magnetic sensor 
mounting position the sensed magnetic field, a noise magnetic field and a signal 
magnetic field as vectors in a coordinate system fixed to the ground. 

The magnetic field generating means for the signal magnetic field in the 
present invention may be a coil. The magnetic field generating means may be 
one electric wire as well, or may also be one electric wire that is straight only in the 
vicinity of the place ol position determination. 

As described above, even if a buried power line, railroad tracks, or similar 
noise magnetic sources are present near a construction site, the present invention 
permits highly reliable position determination without being affected by noise 
magnetic fields generated by such noise magnetic field sources. 

The present invention is intended for measuring (he digging position in the 
non-open-cut method of excavation, but is applicable as well to many technical 
fields that involve position determination by sensing magnetic fields. 
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